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Abstract. We propose a new methodology for probing the cosmological variability of a from pairs of Fell lines 
(SIDAM, single ion differential a measurement) observed in individual exposures from a high resolution spectro- 
graph. By this we avoid the influence of the spectral shifts due to (i) ionization inhomogeneities in the absorbers 
I and (ii) non-zero offsets between different exposures. Applied to the Fell lines of the metal absorption line sys- 

tem at Zabs ~ 1.839 in the spectrum of Q 1101-264 obtained by means of the UV- Visual Echelle Spectrograph 
; (UVES) at the ESQ Very Large Telescope (VLT), SIDAM provides Aa/a= (2.4 ± 3.8stat) x lO"''. The Zate = 

1.15 Fell system toward HE 0515-4414 has been re-analyzed by this method thus obtaining for the combined 
sample Aa/a = (0.4 ± 1.5atat) x 10~^. These values are shifted with respect to the Keck/HIRES mean Aa/a = 
(-5.7± l.l.tat) X lO-'' (Murphy et al. 2004) at very high confidence level (95%). The fundamental photon noise 
QQ , limitation in the A a /a measurement with the VLT/UVES is discussed to figure the prospects for future obser- 

' vations. It is suggested that with a spectrograph of ~ 10 times the UVES dispersion coupled to a 100 m class 

I telescope the present Oklo level (Aa/a > 4.5 x 10~*) can be achieved along cosmological distances with differential 

. measurements of Aa/a . 
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^ 1. Introduction transforms into lAa/aj = |(a^ — a)/a| < 3. 8x 10 ifa^, 

. . the value of a at redshift z, is a linear function of t. The 

> The variability of the fundamental physical constants over functional dependence of the gauge-coupling constants on 

cosmic time, firstly proposed almost seven decades ago ^ -^^ however, unknown and theory predicts even osciUa- 

^ , (Milne 1937; Dirac 1937), has been subsequently mvesti- ^-^^^ ^^^-^^ ^^^^^^ cosmological evolution (e.g., 

; gated m various aspects by many authors (for a review Marciano 1984). In this regard the astronomical observa- 

see, e.g., Uzan 2003). Hons are the only way to test such predictions at different 

The Sommerfeld fine-structure constant, a = /he, space-time coordinates [see, e.g., Mota & Barrow (2004) 

which describes electromagnetic and optical properties of where the effects of inhomogeneous space and time evolu- 

atoms, is the most suitable for time variation tests in both tion of a are studied]. 

laboratory experiments with atomic clocks and astronom- xhe astronomical measurements of the fine-structure 

ical observations. The value of this constant is known with gpUttings of emission lines in distant galaxies started by 

high accuracy, a = 1/137.035 999 76(50) (Mohr & Tailor Savedoff (1956) are summarized in a recent comprehensive 

2000), and its time-dependence is restricted in the labo- ^ork by BahcaU, Steinhardt & Schlegel (2004, hereafter 

ratory experiments^at the level of (M[a{t)]/dt = (-0.9 ± ggS): Aa/a = (0.7 ± 1.4stat) x lO"'^ in the range 0.16 < 

2.9stat) X 10-15 yr-\ corresponding to an upper limit of ^ < 0.80^. The most stringent bound in the overlapping 

\d\n[a{t)]/dt\ < 3.8 x 10'^^ yr'^ (Fischer et al. 2004). At interval z < 0.45, stemming from the radioactive decay 

the cosmological time-scale t - lO^o yr (z > 1), this hmit j-gtes of certain long-lived nuclei found in meteoritic data. 



is set by Olive et al. (2004): Aa/a = (-8 ± 8stat) x 10"^. 

Send offprint requests to: S. A. Levshakov j,^^ y^ie interval comparable with the meteoritic analysis 
lev@astro.ioffe.rssi.ru 



* Based on observations performed at the VLT Kueyen tele- Statistical errors, CTstat = dispersion/y'n, are used in the 

scope (ESO, Paranal, Chile). The data are retrieved from the text to specify the error of the sample mean, if not indicated 
ESO/ST-ECF Science Archive Facility otherwise. 
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Table 1. ESO UVES archive data on the quasar Q 1101-264 
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Cols. 6,7, and 10 list temperature and air pressure inside the spectrograph at the beginning and the end of the exposure 



is covered by the Oklo natural reactor {At ~ 2 x 10^ 
yr). Recent reanalysis of the isotopic abundances in the 
samples taken from Oklo provides an intriguing result that 
the value of a was larger in the past: Aa/a> 4.5 x 10~^ 
(Lamoreaux & Torgerson 2004). 

At higher redshifts the variability of a can be tested by 
observations of small shifts between different ionic tran- 
sitions in the absorption-line spectra of quasars (Bahcall, 
Sargent & Schmidt 1967). This technique, now known as 
the alkali doublet (AD) method, was utilized in numerous 
studies (for a review, see BSS). The best constraint ob- 
tained by the AD method is Aa/a = (—0.5 ± 1.5stat) x 
10-5 (Murphy et al. 2001). 

The AD method was generalized by Webb et al. (1999) 

and Dzuba et al. (1999, 2002) in the many-muhiplet (MM) 
method which provides an order of magnitude improve- 
ment in the accuracy of the estimations of Aa/a. Being 
applied to 143 metal absorption systems (the dominant 
ions Mgii, Fell at z < 1.8, and Alii, Sill at z > 1.8) 
identified in the Keck/HIRES spectra of quasars, the MM 
method indicates, opposite to the recent Oklo result, a de- 
crease of a with cosmic time: Aa/a = (—5.7 ± l.lgtat) x 
10-^ in the redshift range 0.2 < 2; < 4.2 (Murphy et al. 
2004, hereafter MFWDPW). 

A potential concern is that this result includes some 
systematics. Indeed, an increasing accuracy of the Aa/a 
measurements requires a careful consideration of the in- 
trinsic structure of the atomic transitions which is formed 
by the isotope shifts and hyperfine splittings (Levshakov 
1994). For instance, the error of the mean (Ji^/ya/a) = 
1.1 X 10"^ corresponds to the error of the line center 
of ~ 20 m s"^ [see eq.(12) in Levshakov 2004, hereafter 
L04] which is about 40 times smaller than the isotope 
shift between ^^'^'^Mgll transitions 3s 3pi/2,3p3/2, 
Av2i-2Q — 850 m s~^ (DruUinger, Wineland & Bergquist 
1980). Unfortunately, the influence of the isotope shifts 
cannot be well specified since we do not know the isotope 
abundances at different redshifts. If the isotope abundance 



ratio indeed varies with z, the isotope shifts may imitate 
the non-zero Aa/a value (Ashcnfelter, Mathews, & Olive 
2004; Kozlov et al. 2004, hereafter KKBDF). On the other 
hand, at metallicities of Z ~ (0.1 — 1) Zq, typical for the 
QSO systems with low ions, the isotope abundances may 
not differ considerably from terrestrial (Murphy, Webb & 
Flambaum 2003; Chand et al. 2004, hereafter CSPA). 

The influence of unknown isotopic ratio and of another 
source of systematics caused by inhomogencous ionization 
structure within the absorber can be considerably dimin- 
ished if only one heavy element like, e.g.. Fell, is used 
in the Aa/a measurement (L04). In spite of a rather low 
present accuracy of the theoretical calculations of the iso- 
tope shift parameters for atoms with more than one va- 
lence electron (^ 50%, KKBDF), the isotopic effect for 
Fell (seven valence electrons in the configuration 3rf®4s) 
is less pronounced than that for Mg il (one valence electron 
in the configuration 3s) for two reasons: (?) iron is heavier 
and its isotope structure is more compact, and (n) the 
relative abundance of the leading isotope ^^Fe is higher^. 
Therefore systematic shifts in Aa/a due to unknown iso- 
topic compositions should be smaller in the analysis of the 
Fe II data. 

This approach applied to the Fell system identified 
at Zabs= 1-15 in the VLT/UVES spectrum of the bright 
quasar HE 0515-4414 {B = 15.0) gave Aa/a = (-0.4 ± 
1.9stat±2.7sys) X 10~^ (Quast, Reimers & Levshakov 2004, 
hereafter QRL). 

The most stringent limit on the variability of a, 
Aa/a = (—0.6 ± 0.6stat) x 10~^, standard deviation 
(^Aa/a = 4 X 10~^, was claimed by CSPA. Their result 
is based on the VLT/UVES observations of 23 absorption 
systems (0.4 < z < 2.3) toward 18 QSOs. In this study, 
the MM method was applied to a more homogeneous en- 
semble of Mg II, Si II, and Fe 11 lines which are not strongly 



^ The terrestrial isotope ratios are '"'''Fe: '^^Fe: ^''Fe: ^®Fe = 
5.8 : 91.8 : 2.1 : 0.3, and ^''Mg: ^^Mg: ^^Mg = 79 : 10 : 11. 
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Table 2. Atomic data of the Fell transitions", the sensi- 
tivity coefficients Q^, and the isotope mass shift constants 
fcjvis'^- Estimated errors are given in parentheses 
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"based on the compilation of Murphy et al. (2003); 

''defined in Sect. 4; ''calculated by Kozlov et al. (2004); 
''multiplet numbers from Morton (2003) 



2. Construction of the Feii sample 

The most accurate measurements of Aa/a can be carried 
out with unsaturated Fe ii hnes which lie outside the Ly- 
a forest and arc not corrupted by any other absorption 
and/or telluric lines. The components of the Fen lines 
should be clearly distinguished. Besides, a QSO should be 
a bright object to provide a high S/N ratio. These require- 
ments are fulfilled for the 2abs = 1.839 system toward Q 
1101-264 with Zem = 2.145 and V = 16.02, which was dis- 
covered by Osmcr & Smith (1977). The absorption system 
was identified by Carswell et al. (1982) and investigated 
with high spectral resolution by Petitjean, Srianand & 
Lcdoux (2000) and Dcssaugcs-Zavadsky et al. (2003) who 
used the same UVES/VLT data as in the present work. 
This system was also studied with lower resolution in a 
series of works by Boksenberg & Snijdcrs (1981), Young, 
Sargent & Boksenberg (1982), CarsweU et al. (1984), and 
Lanzetta et al. (1987). 



saturated and show a less complex structure as compared 
with the metal profiles from MFWDPW. 

It should be noted, however, that the measurements 
with dispersions Aa/a ^ 2x 10^^ (Aw ^ 60 m ) are al- 
ready at the sensitivity limit of the UVES and, hence, they 
may be affected by the data reduction procedure and/or 
by the changing in time characteristics of the device itself. 
There are three main sources of potential systematic er- 
rors in the absolute velocity scale: {i) temperature and (ii) 
air pressure instability, and (iii) mechanical instabilities 
of unknown origin. For instance, a change of 1 millibar (or 
a change of 0.3°C) induces an error in radial velocities of 
~ 50 m s-i (Kaufer, D'Odorico & Kaper 2004). 

In the present paper we investigate step by step 
the data reduction procedure usually applied to the 
VLT/UVES spectra in order to reveal possible system- 
atics and to design the optimal measurement technique 
to avoid them. The analysis is based on another Fe ii sys- 
tem identified in the absorption-line spectrum of quasar 
Q 1101-264. 

The paper is organized as follows. In Sect. 2 we de- 
scribe observations and data reduction performed to con- 
struct our Fc II sample. Small velocity shifts (equivalent to 
0.1-0.2 pixel size) between the scientific exposures are an- 
alyzed in Sect. 3. Being ignored, such shifts may affect the 
shapes of the Fe ii profiles in the co-added spectra, leading 
to their inconsistency. Sect. 4 presents the key procedure 
used to study the time dependence of a, including error 
estimation. The results of our analysis are explained in 
Sect. 5. Sect. 6 gives the fundamental noise limitation in 
the Aa/a measurement due to photon count and spectral 
profile. The quality factor similar to that introduced by 
Connes (1985) to optimize measurements of the stellar ra- 
dial velocities is defined and calculated in this section. It 
can be used to control the sample dispersion in Aa/a. 
The obtained results and future prospects are discussed 
in Sect. 7. Our conclusions are given in Sect. 8. 



2.1. Observations and data reduction 

The observations were acquired with the UVES at the 
VLT 8.2 m telescope at Paranal, Chile, and the spectral 
data were retrieved from the ESO archive. The high res- 
olution spectra of Q 1101-264 were obtained during the 
UVES Science Verification programme for the study of the 
Ly-a forest (Kim, Cristiani & D'Odorico 2002). The spec- 
tra used here were recorded with a dichroic filter which 
allows to work with the blue and red UVES arms simulta- 
neously as with two independent spectrographs. The stan- 
dard settings with central wavelengths at A437 nm and 
A860 nm were used for the blue and red arms, respec- 
tively. From the blue spectra we used only order 102, and 
from the red spectra - orders 90 and 83, where Fell lines 
suitable for the Aa/a measurement are observed. Details 
of the observations are presented in Table 1. Cols. 1 and 
2 give the exposure archive number and the setting used. 
Cols. 3-7 list the QSO data, and Cols. 8-11 the data of the 
closest in time calibration thorium-argon (TliAr) lamp. 

In this work we analyze five exposures (No. 5/6, 11/12, 
19/20, 21/22, and 23/24) of 3600 s obtained over four 
nights in February 2000 with seeing as indicated in Col. 5 
of Table 1. The slit widths were both set at 0.8 arc- 
sec and the CCDs were read-out in 1x2 binned pixels 
(spatial X dispersion direction). The resulting spectral res- 
olution as measured from the ThAr emission lines is of 
FWHM ~ 6.0 km s'^ in the blue (A ~ 4570 A), and of 
~ 5.4 km s"^ in the red (A ~ 7380 A) 

We used the UVES pipeline (the routines implemented 
in MIDAS-ESO data reduction package for UVES data) to 
perform the bias correction, inter-order background sub- 
traction, flat-fielding, correction of cosmic rays impacts, 
sky subtraction, extraction of the orders in the pixel space 
and wavelength calibration. A modified version of the rou- 
tine 'cchred' of the context ECHELLE inside MIDAS was 
used to calibrate in wavelength the echelle spectra with- 
out rebinning. In this way we have the reduced spectra 
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Fig. 1. Combined absorption-line spectra of Fell associated with the 2;abs= 1.839 damped Ly-a system toward Q 
1101-264 (normalized intensities are shown by dots with Icr error bars). The zero radial velocity is fixed at Zabs = 
1.83888. Smooth lines are the synthetic Voigt profiles convolved with the point-spread spectrograph function. The 
mean signal-to-noise ratio per pixel at the continuum level is indicated in each panel. The normalized Xmin = 3.995 
(the; mimbcr of degrees of freedom ;/ = 272) shows that the Fell profiles obtained in the framework of the standard 
data reduction procedure are not self-consistent 



with their original pixel size in wavelength. This corre- 
sponds to a sampling of 50 mA (3.3 km s^^ /pix) for the 
blue and 55 mA (2.2 km s~^/pix) for the red. Typical 
rms of the wavelength calibration is less than one-fiftieth 
of a pixel, or ffrms ^ 1 mA ( < 60 m s~^ ). The observed 
wavelength scale of each spectrum was transformed into 
vacuum, heliocentric wavelength scale (Edlen 1966). At 
first stage of our analysis we followed the standard proce- 
dure and added together the single extracted spectra using 
weights proportional to the square of their S/N. Before be- 
ing added, the spectra were re-sampled to an equidistant 
wavelength grid (50 mA/pix) using linear interpolation. 

2.2. A test for concordance of the Fe ii profiles 

The Fell profiles in the Zabs = 1.8389 systems toward Q 
1101-264 consist of at least 6 subcomponents spread over 
100 km s~^ (Dessauges-Zavadsky et al. 2003). The central 
component, seen in all but one Fe II lines from Table 2 (the 
undetected A1611 A has the smallest oscillator strength), 
is marginally blended in the blue wing with weak compo- 
nent at u = —17.6 km (see Fig. 1). Henceforth we will 
use this central component in the Aa/a measurement. 

We carefully checked that the profiles shown in Fig. 1 
are free from cosmic rays and telluric absorptions. The 
combined spectra reveal high signal-to-noise ratios, 100 < 
S/N < 160. The continuum levels are very well defined 
for all iron lines since in all cases there are large 'con- 
tinuum windows' at both sides of each Fell absorption. 
But, nevertheless, we failed to find a model which could 
adequately describe the stacked profiles. Using the mul- 



ticomponent Voigt profile fitting procedure we obtained 
the best normalized xt = 3.995 for the number of degrees 
of freedom p = 272. The source of such high xl value is 
clearly seen in Fig. 1: some portions of the Fell profiles 
arc not in concordance with each other. 

This inconsistency in the combined Fell spectra may 
reflect some hidden problems in the standard data reduc- 
tion procedure. We investigate this question in the next 
section. 

3. Velocity shifts between individual exposures 

In order to understand the origin of the discrepancy be- 
tween Fell profiles, we firstly investigated distortion and 
wavelength calibration. It turned out that two Fell lines 
A2344 and A2374 lie close to the starting and end points 
of the echelle orders where distortions in the spectral sen- 
sitivity are the largest. These lines were excluded from the 
further analysis. The remaining Fe II lines are observed in 
the central parts of the echelle orders. 

Then we selected isolated and unblended ThAr lines 
from these central spectral regions and compared their 
profiles. The time intervals between calibration exposures 
are large and the ThAr spectra were taken under difi^erent 
conditions. Namely, the differences between the mean tem- 
peratures and air pressures are AT^/q_ii/i2 = —0.40° C, 
ATn/12-21/22 = 0.15° C, and AP^/e^u/u = -2.73 mb, 
A-Pii/12-21/22 = 1.59 mb, respectively (see Table 1). 
Changes in T and P lead to the velocity shifts between 
exposures. For instance, emission profiles plotted in Fig. 2 
show that in the blue arm AAs-n ~ —6 mA (At; ~ —390 
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m ) and AA21-11 — 4 mA (Ai> ~ 260 m ), whereas 
in the red arm AA22-12 equals AA21-11, as expected, but 
AA6_i2 — —14 rnA (Aw ~ —620 ni ) which is 1.6 times 
larger than ^v^-\\ (exposures 11 and 12 are taken for 
references). If the velocity drifts were entirely due to tem- 
perature and pressure changes, then wc would observe a 
monotonous shift identical in the blue and red arms, as 
in the case of the 21st and 22nd exposures. The deviation 
from this behavior revealed for the 5th and 6th exposures 
indicates that the third factor - some mechanical insta- 
bilities - may also be important. Since it is impossible to 
control what kind of distortions are caused by mechanical 
instabilities, we remove the scientific exposures 5 and 6 
from the study. 

This analysis shows that in order to avoid uncertain- 
ties in the wavelength calibration caused by changes in 
temperature and/or pressure, the calibration ThAr lamps 
must be taken just before and after the scientific expo- 
sure. If for some exposures the relative shifts of the ThAr 
lines in the blue and red arms turn out to be different, 
these exposures must be excluded from the consideration 
since they are influenced by mechanical instabilities which 
cannot be properly corrected. 

Following the standard reduction procedure the re- 
maining scientific exposures should be co-added to en- 
hance S/N. However, this step can introduce additional 
uncertainties because every calibrated exposure has its 
own velocity offset. Note that the combination of the ex- 
posures by means of the MIDAS package requires the same 
starting wavelength, number of points, and the step size. 
Such procedure cancels out the original non-zero offsets 
between the individual exposures and smears out small 
shifts. This can perturb the line centroids, especially those 
of narrow absorption lines with sharp intensity gradients. 
Since in the Aa/a measurement the accuracy of the line 
position is the key factor, it is more appropriate to work 
with individual exposures which have lower S/N rather 
than with a high S/N combined spectrum. 

Accounting for these argTiments wc selected the scien- 
tific exposures with Fell lines shown in Fig. 3. All spec- 
tra are vacuum and heliocentric calibrated with only one 
ThAr spectrum (exposures 11 and 12 for the blue and red 
frames, respectively), they are not combined and not re- 
sampled. The unknown velocity offsets arising from such 
calibration are canceled out in the differential measure- 
ments of Aa/a as described in the next section. 

The key line in our approach is Fell A1608 since its 
sensitivity coeflBcient is negative, while the other iron tran- 
sitions have positive Q- values (Table 2). A relatively low 
strength of the A1608 line requires S/N > 30 in individual 
exposures in order to measure its center accurately. This 
requirement is fulfilled for Q 1101-264 (see Fig. 4, where 
the S/N values are indicated). It is also important that for 
this system the accuracy of the normalization of the Fe II 
profiles leaves no room for doubts since the local continua 
(the horizontal lines in Fig. 3) are calculated with high 
precision. In each panel in Fig. 3, dots with error bars 
indicate the mean intensities and their Icr uncertainties 
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Fig. 2. Emission line profiles from the thorium-argon cal- 
ibration lamp. The lines are placed close to the centers 
of the UVES echellc orders of different settings which 
are indicated at the top left hand corners in the panels. 
The wavelength step size corresponds to the original pixel 
width. The numbering of the histograms arc in accord 
with the list in Table 1. The 11th and 12th exposures are 
taken as reference. Note the difference between the relative 
shifts 5-11 and 6-12, while 21-11 and 22-12 show identical 
velocity offsets 

utilized in the linear regression analysis applied to deter- 
mine the local continuum level. The uncertainty of the 
calculated continuum is found to be less than 1% for all 
selected Fe ii spectra. 

4. The single ion difTerentia! a measurement 

The standard many-multilct (MM) technique (Webb et 
al. 1999, Dzuba et al. 1999, 2002) calculates Aa/a us- 
ing a quite complex procedure. This leads to potential 
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Fig. 3. Unnormalized portions of the Q 1101-264 spectra with the Fell hnes (-Jabs = 1-839) selected for the Aa/a 
measurement. All lines are located near the centers of the UVES echelle orders. The corresponding settings are indicated 
at the bottom left hand corners in each panel. Dots with error bars are the mean intensities and their Icr uncertainties 
- the values used to estimate the local continuum level by means of the linear regression analysis. For the selected 
Fell lines, the continuum level is known with an accuracy better than 1% 



systematic uncertainties which may affect both the line 
profile and the origin of the wavelength scale. The sys- 
tematic effects are thoroughly discussed in Murphy et al. 
(2003), BSS, and L04. In this section we describe how 
to use a modified MM procedure to calculate Aa/a di- 
rectly from the differences between the wavelengths of a 
pair of Fell transitions observed in the individual expo- 
sures. Thereafter this approach is called a 'single ion dif- 
ferential a measurement' (SID AM). An important point 
is that being applied to the lines from the same exposure, 



this method does not depend on the unknown offsets of 

the wavelength scale. 

The MM method utilizes the fact that the energy of 
each line transitions depends individually on a change in 
a. Namely, the frequency of each transition has its own rel- 
ativistic correction to the changes in a which is expressed 
by the coefficient q (Dzuba et al. 1999, 2002). Then, the 
systemic rest wavenumber oJz = 1/Xz '^s given by 



— loq + q (Aa/a) (2 + Aa/a) , 



(1) 
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where = 1/ Aq is the laboratory wavenumber. 

In fact, this approach is similar to the method devel- 
oped by Varshalovich & Levshakov (1993) in order to infer 
the cosmological variability of the proton-electron mass ra- 
tio, ^ — m-p/mc, from the analysis of molecular hydrogen 
H2 absorption lines. The dependence of the frequencies of 
electron- vibro-rotational transitions on a change in /i dif- 
fers for individual transitions and can be characterized by 
the dimensionless sensitivity coefficient /C. The estimation 
of A/x/^ can be obtained from linear regression analysis of 
the {z,]C} pairs from a sample of the H2 lines (Potckhin 
et al. 1998; Levshakov et al. 2002b). Below this technique 
is applied to the Aa/a measurement. 

We can re-write in linear approximation 

(lAa/al < 1) in the form (L04): 



^ Zf) + Ha Qi 



(2) 



where Zi = Aobs,i/Ao,i — 1, Qi ~ (li/^a.i the dimensionless 
sensitivity coefficient, and the slope parameter Kq, is given 

by 



-2(l + ^o)(Aa/a) 



(3) 



If Aa/a is non-zero, Zi and Qi will be correlated. The 
slope Ka and the intercept zq can be found from the linear 
regression analysis of the observed redshifts of the line 
centroids Zi vs. Qi (for more details, see L04 and QRL, 
where this method is called a 'regression' MM analysis, 
RMM). 

In this approach, the value of Aa/a can be directly 
estimated from a pair of lines with different sensitivity 
coefficients. It is easy to show that 



Aa 

a 

or 

Aa 

a 



{Z2 - Zi) 



(2 + Zi+Z2)(Ql- Q2) 



{Z2 - Zi) 



2(l + z)(Qi-Q2) 



(4) 



(5) 



The ratio {z2 — zi)/{l + z) is invariant and does not depend 
on systematic velocity shifts between different scientific 
exposures. 

The computational procedure can be slightly modified 
if we use another variables similar to those introduced in 
the AD analysis by BSS. 

Let us consider an absorption line system at redshift 
z where a set of Fell transitions is formed at the corre- 
sponding cosmic time t. Then, the ratio of the observed 
wavelengths \i{t) to (1 -f z) defines the rest wavelengths 
X'i{t) which may differ from their present-day values Ai(0), 
if Aa/a ^ 0. From one obtains in linear approximation 



A^t) = A,(0)(l-2Q,Aa/a) 



(6) 



R{t) 



Following BSS, we define R{t) and 77 by the relations 
\2{t)-\i{t) 



Ai(i)+A2(i) ' 



(7) 



and 

A2 ^ 

The cosmological redshift of the absorber in the ex- 
pression for R cancels out and thus 



(8) 



A2(i)-Ai(t) _ A^(t)-AUt) 



(9) 



\l{t)+\2{t) \\{t)+X'^{t) 

From 0-©, one obtains after some algebra [taking into 
account ri{t) w ?7(0)] 

Aa 1 



^(^ + 2) A R{t)\ 
a 4(Q2-Qi) ?? + ! V ^(0)/ 



(10) 



In some cases (see Sect. 7.2) this formulation can be more 
convenient for calculations. 



4.1. Error estimations in tine SI DAM technique 

In this section we calculate the uncertainty on individ- 
ual /S.a/a measurement caused by the errors in the wave- 
lengths. The standard method of error propagation is 
used. 

From eq.|01) we deduce (taking into account Z2 ~ zi) 



1 



11) 



where the relative error 5z = (Jz / + z) is equal to 

5.^{5l + 5lS/' . (12) 

If (5ao ^ ^\ and w , then 
1 



"'Act/a 



\/2|Qi- Q2 



(13) 



If the relative errors for both observational and laboratory 
lines are comparable, then the best estimation of Aa/a is 
known with the error 



1 



\Ql-Q2 



(14) 



Using relation H14|) . it is easy to estimate the high- 
est possible precision in the Aa/a measurement with the 
laboratory Fe 11 data at hand. Assuming that the observed 
line centers are known with the accuracy (5a — 5 x 10~* 
(see Table 2), one obtains ctaq/q — 8 x 10~^. This corre- 
sponds to the uncertainty a\ = 0.25 mA (Aw = 15 m s~^ ) 
at 5000 A. The value Aw ~ 15 m s^^ is 4 times lower than 
the limit sensitivity of the UVES, and, hence, with good 
statistics (the sample size > 16) Aa/a can, in principle, 
be probed on the level of < 10^^. 

We note that the best precision reached today in 
the measurements of stellar radial velocities of relatively 
bright stars is a few m s^^ . Observations with HARPS 
provide, for example, the rms uncertainty of 2 m s^^ (Pepe 
et al. 2004) which is close to the HARPS fundamental 
noise limitation. For the Aa/a measurement with UVES, 
the fundamental noise limitation due to spectral profile 
and photon count is considered in Sect. 6. 
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Fig. 4. Individual exposures of the Fell lines in the spectrum of Q 1101-264 (normalized intensities are shown by 
dots with 1(7 error bars) and over-potted synthetic profiles (smooth curves) calculated from the joint analysis. The 
minimization procedure gives x^jj^ = 1.097 per degree of freedom (the number of data points m = 973, the number 
of fitting parameters p = 39). Dashed vertical lines mark positions of the individual components. The mean signal-to- 
noise ratio per pixel at the continuum level is indicated in each panel. Columns from the left to the right correspond 
to the exposures No. 11/12, 19/20, 21/22, and 23/24, in accord with Table 1. Spectral data shown in these columns 
are centered relative to the following fixed redshifts (from the left to the right): 1.8389040, 1.8389041, 1.8389065, and 
1.8389009 



4.2. Errors caused by the isotopic shifts 

The Aa/a measurement on the scale of 8 x 10~^ can be 
affected by the isotopic shifts if the ratio J = ^^Fe/^^Fe 
varies between the absorbers (the influence of the isotope 
^^Fe is negligible because of its relatively low abundance). 
Following KKBDF and using their mass shift constants 



kus for Fell transitions (Table 2), we can estimate the 
shift of the line center of gravity mimicking a non-zero 
Aa/a : 

,„. = ^i„-^-»^.„,(i,-^), (15) 
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Table 3. SID AM analysis: optimized centroid positions of the Fell lines in the main subcomponent and Aa/a 
calculated with ^ or 



Exp. 


Aobs, A 


Aobs, A 


Aobs, A 


Aobs, A 


(Aa/aj 


pair, 


iNO. 


AibUo 


AZooZ 


AzooD 


AzdUU 


units ot iU j 


Al/ A2 


(I) 
('-1 


(2) 


(3) 


(A) 


(5) 


(6) 


(7) 


11/12 


4566.24408 


6764.43913 


7343.27915 


7381.65008 


1.199 


1608/2382 












-2.103 


1608/2586 












0.098 


1608/2600 


19/20 


4566.23755 


6764.44738 


7343.26356 


7381.64359 


-1.047 


1608/2382 












-1.525 


1608/2586 












— U.oyo 


iOUo/ZDUU 


ZL/ ZZ 


4oDD.z4DDi 


D /D4.4oyio 


7343.25994 


7381.63961 


t.DDO 


iOUo/Zooz 












0.539 


1608/2586 












1.860 


1608/2600 


23/24 


4566.25102 


6764.44933 


7343.27203 


7381.65019 


1.209 


1608/2382 












-0.028 


1608/2586 












1.442 


1608/2600 








sample mean (Aa/a): 


0.24 










error of the 


mean a/y/n: 


0.38 





where A' and A are the isotope mass numbers, and x 
(in per cents) shows the reduced fraction of the leading 
isotope A. If, for example, x = 10% (i.e. Iz/Iq — 3), then 
5u)c — —0.004 cm"-'^ for Fell transitions (see Table 1 in 
KKBDF). 

Numerical simulations of the explosive yields show, 
however, that at low metallicities Z < Zq the ratio 
Iz/Iq < 1 (e.g., Chieffi & Limongi 2004). Thus, in the 
high-z absorbers one may expect an approximately con- 
stant isotopic shift SlOc = — 0.12A:msM^ — 0.002 cni^^, 
which is equivalent to the positive shift Aa/a ~ 7.7 x 10~^ 
(note that an offset of the velocity scale of 10-15 m s~^ can 
produce the same effect). However, this shift is canceled 
out in the differential Aa/a measurements of a few high-z 
absorbers having the same metallicities. This can essen- 
tially improve the limiting accuracy of (TAa/a ~ 10~^ set 
by KKBDF for Fe ii samples. 

5. Analysis and the Aa/a results 

In this section we study the Fell profiles selected from 
individual exposures and derive the position of the line 
centroid of the main absorption component seen at zero 
radial velocity in Fig. 4. The total number of the analyzed 
profiles is L = 16. Since the line profiles are very com- 
plex, the distribution of the velocity components along 
the line of sight is crucial for the following Aa/a mea- 
surement. To construct the model for the radial veloc- 
ity distribution, we start with the analysis of the profiles 
from the individual scientific exposures to fix their refer- 
ence frames - the mean Zabs values. The following redshifts 
were determined: Zn/12 = 1.8389040, Z19/20 = 1.8389041, 
Z21/22 = 1.8389065, and Z23/24 = 1.8389009. 

Our model is based on the natural assumption that 
Fell lines have similar profiles, i.e., (1) the number of 



subcomponents Ug is identical for all Fell lines, (2) the 
Doppler bi parameters are identical for the same ith sub- 
components, (3) the relative intensities of the subcompo- 
nents Tij and (4) the relative radial velocity differences 
Avij between the subcomponents are fixed for the ab- 
sorber. We also assume that the main broadening is caused 
by bulk motion. 

Then, the Fell profile is described by the sum of rig 
Voigt functions: 

t(') = A^i ^ r^,l V [{v -V,- Av,,i)/b,\ , (16) 

where Tv is the optical depth at radial velocity v within 
the line t, Ni is the column density of the main compo- 
nent, Ti l — Ni/Ni, V£ is the center of the main component 
in the line £, and Avi^i — 0. 

Our model is fully defined by specifying A^i, 
{Awi_i}"^7\ and {vi}^^-^. All these param- 

eters are components of the parameter vector 6 = 
{61, 6*2, ... , Op}- To estimate 9 from the Fell profiles, we 
minimize the objective function 

L nil 

^'(^) = ^ E E l^^y^ - -^"'f /^L- ' (17) 

where T^^^ is the observed normalized intensity of the 
spectral line £, aej is the experimental error within the 

jth pixel of the line profile, and v = ™^ ~ pj 

number of degrees of freedom. J-'^^j{9) is the calculated 
intensity convolved with the spectrograph point-spread 
function. 

The synthetic profiles for the best are shown by 
the smooth lines in Fig. 4. The optimal number of sub- 
components is ris = 8. Their positions are marked by 
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Fig. 5. Upper panel: The residuals e = (JF^^' — J^°^^^ for 
all Fell profiles from Fig. 4. Dashed vertical lines indicate 
the positions of the subcomponents. Lower panel: Shown 
are the distribution of e (histogram) and the overplotted 
normal distribution (smooth curve) with the sample mean 
(e) = 3.5 X 10~^ and the dispersion Ue = 2.9 x 10~^ esti- 
mated from the total data set except for two 'hot pixels' 
with e < —0.1 (seen in the inset) 



vertical dotted lines in the panels in Fig. 4. The corre- 
sponding Xmin ~ 1-21 (for m = 975 and p — 39) is, 
however, too high: with i/ = 936, the expected mean is 
xt = 1±0.05 (ItT c.L). However, the analysis of the resid- 
uals e — (JF'^^' — jF°t>s-j sj^own in Fig. 5 reveals two 'hot 
pixels' with e < —0.1 which deteriorate this Xmin value'^. 
After removing these points, we find Xmin ~ 1-097 (for 
m = 973 andp = 39). 



^ The first 'hot pixel' with e — —0.15 is seen in Fig. 4, expo- 
sures 23/24, panel Fell A2586.6, at ?; = —66.6 km s~^ , whereas 
the second one with e — —0.24, the same exposures, panel Fen 
A1608.5, V = 65.6 km s~^ is beyond the frame. 



The estimated wavelengths of the main Fell compo- 
nents are listed in Table 3. These are the best fitting quan- 
tities. We do not calculate their errors since the spectral 
data from different exposures show very similar S/N ra- 
tios which allow us to calculate the sample mean (Aa/a) 
without weights. The results are given in Col. 6 of Table 3. 
Col. 7 lists the Feii pairs used in a particular Aa/a mea- 
surement. We find (Aa/a) = (2.4 ± 3.8) x lO'^, and the 
corresponding root mean square CTrms = 1.3 x 10^^. 

6. The fundamental noise limitation in the Aa/a 
measurement 

In Sect. 4, we described the computational procedure used 
in this study - the single ion differential a measurement, 
SIDAM. A simple form of the deduced eqs. Q and IpU 
allows us to compute directly the fundamental uncertainty 
in the Aa/a measurement due to photon noise. This anal- 
ysis is based on the results obtained by Connes (1985) and 
by Bouchy, Pepe & Queloz (2001) who calculated the fun- 
damental noise limitation in the Doppler shift measure- 
ments. 

Let us consider a digitalizcd and calibrated spectra of 
a pair of Fe ii lines which are obtained with a high stabil- 
ity spectrograph. Let AApix be the pixel size (the wave- 
length interval between pixels). Assume further that the 
spectrograph point-spread function can be described by 
a Gaussian with FWHMgp = 2 AApix (the Nyquist limit). 
The observed Feii lines are supposed to be isolated, sin- 
gle component, and unsaturated. Their apparent width, 
FWHMiinc, is caused by the convolution of the spectro- 
graph point-spread function with the 'true' profile. The 
width of the true profile is defined by the quadratic sum 
of the thermal and turbulent components: FWHM^j.^^, = 
FWHM^j^^^^ + FWHM^^^^. Since Fell hues are usually 
observed in the damped Ly-a systems where kinetic tem- 
perature is low (~ 100 K, FWHMthcrm ^ 0.3 km s'^ ) 
and the turbulent broadening is a few km s^^ , FWHMtruo 
might be less or about FWHMgp. 

The error in the line center caused by counting statis- 
tics is given by Bohhn et al. [1983, eq.(A15)]: 



AA, 



(18) 



where Wobs is the observed equivalent width, J\fc is the 
mean number of photoelectrons per pixel at the continuum 
level, and M is the number of pixels covering the line 
profile^. 

For Gaussian profiles, M can be equal to 
2.5FWHMiinc/AApix, producing 

A/ = 5e , (19) 

where C = [1 + (FWHMtrue/FWHMsp)^ ]i/2. 

^ The second term in Bohlin's et al. eq.(A15) is neglected 
since we assume the lines are rather strong (the apparent cen- 
tral optical depth ro ~ 1), and isolated, i.e. the local continuum 
level is known with a sufficiently high precision. 
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Table 4. Fundamental photon noise limits cr^^j^ (in units 
of 10~^) for Fell pairs of lines (A1608/X) from the Zabs = 
1.839 system toward Q 1101-264. The exposure time icxp 
and the total efficiency etot are set to 3600 s and 0.15, 
respectively 



Line: 


A1608 


A2382 


A2586 


A2600 


Wobs, A 


0.0446 


0.2779 


0.1239 


0.2770 


AApix, A 


0.05 


0.06 


0.05 


0.05 




1.83 


1.83 


1.93 


1.93 


Q-factor 


0.11 


0.58 


0.29 


0.64 


a^, mA 


5.7 


1.3 


2.3 


1.0 


-lim 




1.1 


1.0 


1.1 



By analogy with the Connes procedure, we can char- 
acterize a line profile by a dimensionless quality factor Q 
which is independent on the flux: 



Q = 



bs 



AApix 5%/5e\/? ' 
and re-write (|18|1 in the form 

AAoix 



Q^/K 



(20) 



(21) 



This equation shows that the error ax decreases almost 
quadratically with decreasing wavelength bin per pixel 
(W^obs does not depend on the spectral resolution, and 
^ - V2, if FWHMsp ~ FWHMtruc) under the assumption 
that Nc is keeping constant. The quality factor for Fell 
lines from the Zabs = 1.839 system ranges between 0.1 and 
0.7 (Table 4). 

The total number of photoelectrons can be estimated 
from the specific fiux — 3.81 x 10"^" ergs s~^ cm~^ 
Hz^^ of a mv = star outside the Earth's atmosphere^: 



by: 



1.05 X 10"^ photon cm^ -2 A-i. ThenWo IS given 



FY Stc\ £tot ^cxp AApix 
2Q0.4 mv 



(22) 



Here stei is the telescope area in cm^; etot the overall detec- 
tion efficiency of the telescope, spectrograph and detector, 
corrected for the contribution of the atmosphere; texp the 
exposure time in s; mv the visual magnitude of the quasar. 

At 5500 A, the UVES efliciency etot ^ 0.15 (Kaufer et 
al. 2004), and with a 3600 s exposure [Nc — 5670), one 
expects a fundamental uncertainty of about 1-5 mA in ax 
for a mv = 16 QSO. Then, fTT|l provides, respectively, the 
fundamental noise of ~ (0.2 — 1.0) x 10~^ in (TAa/a for 
one measurement of a Fell pair of lines. 

The fundamental photon noise limits for Fell pairs 
from the Zabs = 1.839 system are given in Table 4. The 
photon noise is computed with the best-fitting parame- 
ters iVi = 1.335 X 10" cm~2 ^nd bi = 5.25 km s"! . 



7. Discussion and Outlook 

7.1. /S.a/a from different studies 

Measurements from the present paper can be combined 
with the previous Fe ii sample from the Zabs =1.15 system 
toward HE 0515-4414 (QRL) to increase statistics. The 
normalized distribution (/•dAa/a= 1) of the resulting 
35 Aa/a values is plotted in Fig. 6 (histogram) along with 
two other recently published results of MFWDPW and 
CSPA which are shown by the dashed and dotted curves, 
respectively, assuming that the measured Aa/aare nor- 
mally distributed with the sample means and standard 
deviations published in these papers. The vertical lines in 
this figure mark the centers of the corresponding distri- 
butions. The distribution shown by the histogram has the 
sample mean (Aa/a) = (0.4 ± 1.5stat) x 10~^ and the 
median (Aa/Q;)incd = 0.1 x 10~^. 

It is to be noted that at a given redshift the sam- 
ple mean (Aa/a) should be the same within the un- 
certainty interval independently on the method or the 
sample used - provided the data are free from any sys- 
tematics. The results presented in Fig. 6 show, however, 
that (Aa/a)Kock/HiRES 7^ (Aa/a)vLT/uvES- The sam- 
ple means of CSPA and our Fell ensemble are in good 
agreement but they differ from that of MFWDPW at the 
95% significance level according to the t-test®. This dis- 
crepancy points to the systematic shift which mimics the 
effect of varying a in the Keck/HIRES spectra. To clarify 
the origin of this systematic shift, one needs more accu- 
rate measurements, which can be carried out with higher 
spectral resolution and with more homogeneous samples. 

The comparison of the distribution widths in Fig. 6 
reveals that the standard deviation in the CSPA sample 
is exceptionally small. For example. Fig. 1(b) in CSPA, 
where the accuracy of wavelength calibration is checked 
through the relative velocity shifts, Au, between the Fe ii 
A2344 and A2600 lines^, shows the dispersion of CAt, — 
0.4 km s^^ . This uncertainty in wavelength calibration 
transforms into the error (Taq/q ^ 2 x 10^^ [see eq.(12) in 
L04] , i.e. , in order to reach the error of the mean cr(Aa/a) ^ 
0.6 X 10~^ (CSPA), one needs a sample of the size n ~ 
1100, which is not the case. Thus, the error of the mean 
(T^Aa/a) estimated by CSPA is in some disagreement with 
their Fig. 1(b). The scatter of Aa/ain the Keck sample is 
about 2 times the CTrms value of our combined Fe ii sample. 

7.2. Future prospects 

Future measurements of Aa/afrom astronomical obser- 
vations should be carried out with high precision to check 
the present Oklo result that a was larger in the past. 



For B and R filters, the corresponding intensities are 
1.52 X 10^ and = 0.65 x 10^ photon cm'^ s'^ A'^ 



The application of the f-test to two normally distributed 
and independent random variables with different variances is 
described in Bol'shev & Smirnov (1983). 

^ The mean (Au) should be consistent with zero in the case 
of good calibration since the sensitivity coefficients Q(Fell 
A2344) ~ Q(Feii A2600). 
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Fig. 6. The distribution of Aa/a (histogram) from the 
^abs = 1.839 and Zabs = 1-15 Fell systems toward Q 1101- 
264 and HE 0515-4414, respectively. The distributions of 
Aa/a from MFWDPW (dashed curve) and CSPA (dot- 
ted curve) are overplotted to demonstrate the discrepancy 
between the Keck/HIRES and VLT/UVES sample mean 
(Aa/a ) values. Both MFWDPW and CSPA distributions 
are assumed to be normal with the means and variances 
taken from the corresponding papers. The vertical lines 
mark the centers of the distributions. The revealed Keck- 
VLT discrepancy of (Aa/a) is statistically significant at 
the 95% confidence level according to the i-test 

Aa/a> 4.5 x 10^* (Lamoreaux & Torgerson 2004). To 
reach the Oklo scale of a- variation, one needs the accuracy 
of at least (TAa/a ~ 10~^, which is 8 times higher than the 
accuracy set by the errors in the laboratory wavelengths 
(Sect. 4.1). At first glance the Oklo level is unachiev- 
able in astronomical observations. However, applying the 
SID AM method to a few Fell systems observed at dif- 
ferent redshifts, one can omit the laboratory wavelengths 
from Aa/a calculations. Such self- calibrating procedure, 
described by BSS, implies that the measured values of 
R{t), eq.JTI), can be fitted to a linear function of cosmic 
time given in the form [cf. eq.(6) in BSS]: 

R{t) = R{0) (1 + StHo) , (23) 

where R{0) is any convenient constant, Hq is the Hubble 
constant, and the slope S is equal in our case to 

' da\ 

H) ■ 



5* = 



(24) 



The constant n is defined through ea. H10() : k — 4(Qi 
Q2)(?? + l)/??('? + 2). For the pair Fell AA1608,2586, n 
-1/4. 



Equation H21|l shows that for narrow lines (^ < ^/2) the 
error (j\ decreases almost quadratically with decreasing 
wavelength bin per pixel if Nc is fixed. Therefore, if an 
efficient spectrograph with ten times the UVES dispersion 
and superior stability can be coupled to a 100 m class 
telescope, one would expect a^™ '-^ 0.03-0.05 mA, and, 
thus, the precision of ~ 10~^ in the Aa/a measurements 
can be achieved. Then, with good statistics the Oklo result 
can be checked at different redshifts. 

If a were indeed larger in the past, then the logarithmic 
derivative of a{t) in ea. (|24|l is negative (since t decreases 
with increasing redshift) and we would observe a positive 
slope S. Otherwise, if a were smaller in the past, S would 
be negative. 

8. Conclusions 

The main results of the present paper are as follows: 



1, 



A data reduction procedure is designed to control the 
accuracy of the wavelength scale calibration and to 
bypass the influence of the spectral shifts due to in- 
stabilities in the instrument which perturb the narrow 
absorption-line profiles. 

A single ion method is proposed to minimize the in- 
fluence of the systematic shifts of the line centroids to 
the accuracy of the Aa/a measurements. 
The mean value (Aa/a) = (2.4± 3.8stat) x 10"^ is ob- 
tained from the analysis of the Zabs = 1.839 Fe ii system 
toward Q 1101-264. Combination of this measurement 
with the Zabs — 1.15 Fe ii system toward HE 0515-4414 
(QRL) increases the sample size to 35 Fell pairs and 
provides (Aa/a) = (0.4 ± 1.5stat) x lO'^. 
We confirm that the VLT/UVES estimations of 
Aa/a differ significantly from the Keck/HIRES result, 
and that the systematic shift is present in the Keck 
data with a probability of 95%. 

The fundamental noise limitation in the 
Aa /a measurement is defined and calculated for 
the Zabs = 1.839 Fell system. It is shown that 
a typical VLT/UVES fundamental noise limit is 
1 X 10"^ for a pair of Fell lines observed 



lim 



in a spectrum of a relatively bright QSO {V — 16.0). 
This estimation is robust since the inclusion of other 
metal transitions (like Mgii AA2796, 2803; Alii A1670; 
Sin AA1526, 1808) cannot improve significantly the 
due to lower sensitivity to the varia- 



error a 



lim 

Aa /a. 



tion of a. From this point of view the dispersion 
"■Aa/a ~ 0.4 X 10^^ and the error of the mean 
""(Aq/q) = 0-6 X 10~^ found in CSPA seems to be 
several times underestimated. 
6. We suggest that with a spectrograph of ^ 10 times the 
UVES dispersion and superior stability the accuracy 
high enough to probe the Oklo result, Aa/a > 4.5 x 
10~^, can be attained in future at new giant telescopes. 
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